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Abstract. The "P32" Astronomical Observation Template (AOT) provided a means to map large areas of sky 
(up to 45x45 arcmin) in the far- infrared (FIR) at high redundancy and with sampling close to the Nyquist 
limit using the ISOPHOT C100 (3x3) and C200 (2x2) detector arrays on board the Infrared Space Observatory 
(ISO). However, the transient response behaviour of the Ga:Ge detectors, if uncorrected, can lead to severe 
systematic photometric errors and distortions of source morphology on maps. We describe the basic concepts of 
an algorithm which can successfully correct for transient response artifacts in P32 observations. Examples are 
given to demonstrate the photometric and imaging performance of ISOPHOT P32 observations of point and 
extended sources corrected using the algorithm. For extended sources we give the integrated flux densities of the 
nearby galaxies NGC 6946, M 51 and M 101. and an image of M 101 at 100 fim. 
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1. Introduction 

From the point of view of signal processing and photom- 
etry, diffraction-limited mapping in the FIR with cryo- 
genic space observatories equipped with photoconductor 
detectors poses a particular challenge. In this wavelength 
regime the number of pixels in detector arrays is limited in 
comparison with that in mid- and near-IR detectors. This 
means that more repointings are needed to map structures 
spanning a given number of resolution elements. Due to 
the logistical constraints imposed by the limited opera- 
tional lifetime of a cryogenic mission, this inevitably leads 
to the problem that the time scale for modulation of il- 
lumination on the detector pixels becomes smaller than 
the characteristic transient response timescale of the de- 
tectors to steps in illumination. The latter timescale can 
reach minutes for the low levels of illumination encoun- 
tered on board a space observatory. The result is that the 
signals from the detectors can depend as much on the il- 
lumination history as on the instantaneous illumination. 
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Unless corrected for, the transient response behaviour of 
the detectors lead to distortions in images, as well as to 
systematic errors in the photometry, especially for discrete 
sources appearing on the maps. In general, these artifacts 
become more long lived and more difficult to correct for at 
fainter levels of illumination, since the transient response 
timescales increase with decreasing illumination. 

For mapping observations in the FIR, the Infrared 
Space Observatory (ISO; Kessler et al. 1996) was equipped 
with two Gallium-doped Germanium photoconductor de- 
tectors which comprised part of the ISOPHOT instrument 
(Lemke et al. 1996). The ISOPHOT C200 detector was a 
2x2 (stressed) pixel array operating in the 100 — 200 /mi 
range, and the ISOPHOT C100 detector was a 3 x 3 pixel 
array operating in the 50 — 100 /im range. Compared to 
the IRAS survey detectors, the ISOPHOT-C detectors had 
relatively small pixels designed to provide near diffraction 
limiting imaging. ISOPHOT thus generally encountered 
larger contrasts in illumination between source and back- 
ground than IRAS did, making the artifacts from the tran- 
sient response more prominent, particularly for fields with 
bright compact sources and faint backgrounds. 

A further difficulty specific to mapping in the FIR with 
ISO was that, unlike IRAS, the satellite had no possi- 
bility to cover a target field in a controlled raster slew 
mode. This limited the field size that could be mapped 
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using the spacecraft raster pointing mode alone, since the 
minimum time interval between the satellite fine point- 
ings used in this mode was around 8 s. This often greatly 
exceeded the nominal exposure time needed to reach a 
required level of sensitivity (or even for many fields the 
confusion limit). Furthermore, the angular sampling and 
redundancy achievable using the fine pointing mode in 
the available time was often quite limited, so that com- 
promises sometimes had to be made to adequately extend 
the map onto the background. 

A specific operational mode for ISO - the "P32" 
Astronomical Observation Template (AOT) - was devel- 
oped for the ISOPHOT instrument to alleviate these ef- 
fects (Heinrichsen et al. 1997). This mode employed a 
combination of standard spacecraft repointings and rapid 
oversampled scans using the focal plane chopper. The 
technique could achieve a Nyquist sampling on map areas 
of sky ranging up to 45x45 arcmin in extent (ca. 70x70 
FWHM resolution elements) on timescales of no more 
than a few hours. In addition to mapping large sources, 
the P32 AOT was extensively used to observe very faint 
compact sources where the improved sky sampling and re- 
dundancy alleviated the effects of confusion and glitching. 

In all, over 6% of the observing time of ISO was de- 
voted to P32 observations during the 1995-1998 mission, 
but the mode could not until now be fully exploited sci- 
entifically due to the lack of a means of correcting for the 
complex non-linear response behaviour of the Ge:Ga de- 
tectors. Here we describe the basic concept of an algorithm 
which can successfully correct for the transient response 
artifacts in P32 observations. This algorithm forms the 
kernel of the "P32Tools" package, which is now publically 
available as part of the ISOPHOT Interactive Analysis 
package PIA (Gabriel et. al 1997; Gabriel & Acosta-Pulido 
1999). Information on the algorithm, as well as the first 
scientific applications, can also be found in Tuffs et al. 
(2002a,b) and Tuffs & Gabriel (2002). The user inter- 
face of P32Tools that connects the P32 algorithm to PIA 
is described by Lu et al. (2002). After a brief overview 
of relevant aspects of the P32 AOT in Sect. 2, we de- 
scribe the semi-empirical model used to reproduce the 
transient response behaviour of the PHT-C detectors in 
Sect. 3. Sect. 4 describes the algorithms used to correct 
data. Examples demonstrating the photometric and imag- 
ing performance of ISOPHOT P32 observations are given 
in Sect. 5, based on maps corrected using the P32 algo- 
rithm. 

2. The P32 AOT 

The basic concept of the P32 AOT, in which imaging was 
achieved using a combination of standard spacecraft re- 
pointings and rapid oversampled scans using ISOPHOT's 
focal plane chopper, is summarised in Fig. 1. There were 
two basic requirements, outlined in the original proposal 
for the AOT (Tuffs & Chini 1990) which led to the final de- 
sign: The first goal was to give ISOPHOT the capability of 
achieving a Nyquist sky sampling (AO = A/2D = 17arcsec 
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Fig. 1. A schematic representation of the P32 observing 
technique. The sky is sampled using a combination of a 
coarse spacecraft raster in the satellite Y and Z coordi- 
nates (horizontal and vertical directions on the figure) and 
finely sampled sweeps of the focal plane chopper in the Y 
coordinate. The spacecraft raster (in this example a 4 x 3 
raster) starts at the maximum value of Z and the min- 
imum value of Y and proceeds as shown by the arrows. 
The position angle (3 of the first scan leg in Y according 
to normal astronomical convention (taken positive E from 
N) is related to the satellite roll angle a by (3 = a + 90°. 
At each spacecraft pointing direction a chopper sweep of 
length 180 arcsec (the diameter of the unvignetted field 
of view of ISO) is made with a sampling of 1/3 of the 
detector pixel separation. This results in chopper sweeps 
of 13 positions separated by 15 arcsec for observations us- 
ing the C100 detector array (as depicted in this example), 
and chopper sweeps of 7 positions sampled at 31 arcsec 
for the C200 detector array. The spacecraft raster sam- 
pling interval in Y is set to be a multiple of the chopper 
sampling interval (in the example depicted it is set to 6 
chopper sampling intervals). A series of uniformly sam- 
pled linear scans for each detector pixel in Y is obtained 
in this fashion, built up from a sequence of registered and 
overlapping chopper sweeps obtained on each spacecraft 
scan leg in Y. The sampling of the target field in Z is con- 
trolled entirely by the spacecraft raster sampling interval 
in Z. In the example depicted the C100 detector detector 
array is stepped in intervals of 1.5 times the detector pixel 
separation in Z. For this particular case, the linear scans 
from different CI 00 detector pixels can be combined to 
sample the sky at half the detector pixel separation in Z, 
resulting in an overall sampling of 15x23 arcsec over the 
target field. 



at 100 /*m) on extended sources subtending up to 50 res- 
olution elements in a tractable observation time. As well 
as providing a unique representation of any arbitrary sky 
brightnesss distribution, this also improved the effective 
confusion limit in deep observations, which, in directions 
of low background was ~ 100 and 20mJy rms at 200 and 
100 ^m, respectively, for ISOPHOT-C. 
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The second requirement was to achieve a redundancy 
in the observed data. This led to the sampling of each 
of a given set of sky directions by each detector pixel at 
different times. The redundancy was achieved in part by 
making several (at least 4) chopper sweeps at each space- 
craft pointing direction. Further redundancy was achieved 
by having an overlap in sky coverage of at least half the 
chopper sweep amplitude between chopper sweeps made 
at successive spacecraft pointing directions in the Y space- 
craft coordinate. Thus, a given sky direction was sampled 
on three different time scales by each detector pixel: on 
intervals of the detector non-destructive read interval, on 
intervals of the chopper sweep, and on intervals of the 
spacecraft raster pointings. Although the redundancy re- 
quirement was originally made to help follow the long 
term trends in detector responsivity during an observa- 
tion, in practice it also proved very useful in removing ef- 
fects caused by cosmic particle impacts (Gabriel & Acosta- 
Pulido 2000) , the so-called "glitches" , from the data. This 
was especially important to optimise the sensitivity of 
deep observations using the CI 00 detector, which are gen- 
erally limited by glitching rather than by confusion. 

Both these requirements led to dwell times per chop- 
per plateau which could be as low as 0.1 s. For a given sky 
brightness, this also had the consequence that the detector 
was generally read faster than in the other AOTs using the 
PHT-C detectors, which allows a higher time resolution in 
the deglitching procedures. A further beneficial side effect 
of the P32 AOT was that because the transient response 
timescale of the detectors is longer than the dwell time of 
the chopper, even for bright sources, the effective limiting 
source brightness that could be observed without saturat- 
ing the detector was higher than for the other AOTs. 

2.1. The P32 "natural grid" 

The fact that the spacecraft pointing increment in Y was 
constrained to be a multiple of the chopper step interval 
means that a rectangular grid of directions on the sky can 
be defined such that all data taken during P32 measure- 
ment fine pointings will, to within the precision of the 
fine pointings, exactly fall onto the directions of the grid. 
This so-called P32 "natural grid" allows maps to be made 
without any gridding function, thus optimising angular 
resolution. The P32 natural grid is also a basic building 
block of the algorithm to correct for the transient response 
behaviour of the detector, as the algorithm solves for the 
intrinsic sky brightness in MxN independent variables for 
a grid dimension of MxN. Here M and N are the number 
of columns and rows on the P32 natural grid. An example 
of the almost perfect registration of data onto the grid is 
shown in Fig. 2. 

All the example maps shown in this paper will be sam- 
pled without a gridding function on the P32 natural grid. 
This means that the data in each of the map pixels is 
independent. This is unlike the situation for most other 
maps calculated using PIA, where a gridding function is 



employed and so raw data can contribute to more than 
one map pixel. 
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Fig. 2. Pointing directions observed towards Ceres at 
105 fim. The spacecraft raster dimension was 3 x 3 in 
the spacecraft coordinates Y x Z. Sky sampling in Y is 
achieved through a combination of chopper scans in Y 
and repointings of the spacecraft. In Z the sky sampling 
is determined alone by the spacecraft repointing interval. 
Each pointing direction for 11833 individual data sam- 
ples for the central pixel of the CI 00 detector array is 
plotted as a cross. The rectangular grid is the "P32 natu- 
ral grid" (see text) for this observation, on which the sky 
brightness distribution is to be solved. The grid sampling 
is 14.99 x 45.98 arcsec. Only points with the on target flag 
set (i.e. not including slews) are plotted. The pointings 
typically lie within 1 arcsec of the centre of each pixel of 
the P32 natural grid. 



3. Transient response of the PHT-C detectors 

Gallium doped germanium photoconductor detectors ex- 
hibit a number of effects under low background condi- 
tions which create severe calibration uncertainties. The 
ISOPHOT C200 (Ge:Ga stressed) and especially the C100 
(Ge:Ga) detectors have a complex non-linear response 
as a function of illumination history on timescales of 
~ 0.1 — 100 sec, which depends on the absolute illumina- 
tion as well as the changes in illumination. The following 
behaviour is seen in response to illumination steps (see 
also Fouks 1992, Acosta-Pulido et al. 2000 and Laureijs et 
al. 2000, Sect 4.2.2.): 

— a rapid jump in signal, which however undershoots the 
asymptotic equilibrium response (often severely for the 
C100 array) 

— a hook response (overshooting or undershooting on in- 
termediate timescales of ~ 1 — 10 s) 

— a slow convergence to the asymptotic equilibrium re- 
sponse of the detector on timescales of ~ 10 — 100 s 

An example showing these effects for a staring observa- 
tion is given in Fig. 3. Also seen are points discarded by 
the deglitching algorithm in PIA, marked by stars. Many 



Tuffs & Gabriel 2003: Photometric Mapping with ISOPHOT 



U 2.2 - 
> 






2.00 
1.90 
1.80 
1.70 
1.60 




1 


1 . 


, i , 


2 4 6 

, 1 


8 10 12 14 

1 , 


60 


80 


100 


120 



Time [si 



Fig. 3. Example of the transient response of a Ge:Ga de- 
tector following an upwards step to a constant level of 
illumination at 100 /im. The hook response is shown in 
the inset. 



glitches are accompanied by a longer lived "tail". It is 
these tails that determine the ultimate sensitivity of the 
ISOPHOT C100 array. 

3.1. Semi-empirical detector model 

Unlike the ISOCAM and ISOPHOT-S Si:Ge detectors on 
board ISO, which operate in the mid-IR, the transient re- 
sponse behaviour of the Ge:Ga detectors of ISOPHOT-C 
are not described by the analytical Fouks Schubert model 
(Fouks & Schubert 1995; Schubert et al. 1995; Acosta- 
Pulido 1998; see Coulais & Abergel 2000 for its applica- 
tion to ISOCAM). Although numerical simulations of the 
transient behaviour of Ge:Ga detectors have been made 
(Hagel et al. 1996), these do not cover the wide range of 
operating conditions experienced by the ISOPHOT detec- 
tors. However, the effects observed in ISOPHOT-C can to 
some extent be qualitatively understood in terms of the 
build up of charge on the contacts between semiconduc- 
tor and readout circuit, as analytically modelled by Sclar 
(1984). Sclar's model can be formulated as the superposi- 
tion of three exponentials to describe the response of the 
detector to a step in illumination. A superposition of three 
exponentials has also been proposed by Fujiwara et al. 
(1995) and Church et al. (1996) to represent the transient 
behaviour. The three timescales involved are a minimum 
requirement to reproduce the hook response and the long 
term approach to the asymptotic response. However, no 
theory exists to predict values for the constants in this 
representation, and thus they must be found empirically. 

In practice it is a formidable and unsolved problem 
to extract the constants from in flight data for as many 
as three temporal components to the transient response 
behaviour. In this paper we describe an approximate so- 
lution involving a superposition of two exponentials, each 
with illumination-dependent constants. This was found to 
adequately represent the detector response on timescales 



greater than a few seconds, though it only approximately 
model the hook response (see also Blomme & Runacres 
2002). Our model can be viewed as a generalisation of the 
single exponential model developed for correction of longer 
stares by Acosta-Pulido, Gabriel & Castaheda (2000). 
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Fig. 4. A schematic representation of the response of the 
semi-empirical detector model used in the P32 algorithm 
to an upwards step in illumination from Soop to Soo- The 
overall response S is the sum of a "slow" signal response 
Si and a "fast" signal response 5*2 (see Sect. 3.1). Si p and 
S 2p respectively denote the signal response of these two 
components immediately prior to the step in illumination. 



In the model considered here, schematically illustrated 
in Fig. 4, the basic detector response S to an instantaneous 
step in illumination from Soo P to Soo is given by the sum 
of a "slow" signal response Si and a "fast" signal response 
S 2 as: 

S = Si + S 2 (1) 

51 = (1 - (3 2 ) Soo (1 ~ exp[-t/n]) + Soi exp[-t/ n ] (2) 

5 2 = ft Soo (1 - exp[-t/r 2 ]) + S 02 exp[-t/r 2 ] (3) 

Soi and So 2 are signals corresponding to the slow and fast 
response components immediately following the illumina- 
tion step. Soi and So 2 are related to the corresponding sig- 
nals immediately prior to the illumination step, Si p and 
S 2p , which contain the information about the illumination 
history. Soi and So 2 are specified by the jump conditions: 



501 — Pi (Soo — Soop) + Sip 

50 2 = S 2p 



(4) 
(5) 



The model is thus specified by four primary parameters: 
the slow and fast timescales t\ and t 2 , the "jump" factor 
Pi and a constant of proportionality for the fast response 
/3 2 . In order to obtain satisfactory model fits to complex 
illumination histories, each of these primary parameters 
had to be specified as monotonic functions of illumination, 
each with three subsiduary parameters: 

Pi = fto + 0u * (6) 
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Table 2. Model parameters for the C200 detector 
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t\ = no + Tn * S^ 7 " 12 

ft = /?20 + 021 * 

T2 = T 2 + T21 * S 1 " 1 " 22 



(7) 
(8) 
(9) 



Thus 12 subsiduary parameters in all are needed to char- 
acterise the transient response behaviour of the detectors. 
In addition, the initial starting state of the detector, given 
by the values of S\ p and S^p immediately prior to the start 
of the mapping observation, must be specified. 
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Fig. 5. Example of a model fit (curve) to the response of 
the central pixel of the C100 array to a constant FCS illu- 
mination starting at t = s. Prior to the FCS illumination 
the detector was viewing blank sky. The horizontal line in- 
dicates the fitted value of the illumination. No data was 
recorded in the first two seconds following FCS switch-on. 



3.2. Determination of model parameters 

The parameterisation of the illumination dependence of 
the primary parameters n, (3\, t 2 and /3 2 was made in en- 
gineering units (V/s), rather than in astrophysical units 
such as MJy/sr. Thus, it was assumed that there is no 



dependence of the detector response on the wavelength of 
the IR photons reaching the detector. Since the individual 
pixels of the C100 and the C200 arrays behave as indepen- 
dent detectors, each of the 12 parameters of the detector 
model had to be separately determined for each detector 
pixel. 

The "slow" primary parameters n and (5\ were found 
from fits to the detector response to fine calibration source 
(FCS) exposures. To make this determination independent 
of the determination of the "fast" parameters, only data 
points corresponding to times of 5s or longer after the 
switch-on of the FCS was used. An example of a fit to the 
signal from an FCS exposure is given in Fig. 5. The model 
fit to the signal response on time scales ranging from a few 
seconds up to the duration of the FCS exposure depends 
mainly on the values of n and /3i . 

A database of values for n and (3\ was built up for a 
sample of FCS measurements, selected so as to correspond 
to the complete range of sky and source illuminations en- 
countered by ISOPHOT during the ISO mission. These 
allowed the dependence of n and [3\ on illumination to be 
determined according to Eqns. 6 & 7. An example showing 
the illumination dependence of (3\ for the central pixel of 
the C100 array is given in Fig. 6a. The solid line is the fit- 
ted illumination dependence of [3\ used to find the values 
of the parameters /3io, /?n and /3i 2 used in the model. A 
corresponding plot for the illumination dependence of n 
is given in Fig. 6b. 

At shorter timescales the response is primarily de- 
termined by the "fast" primary parameters t 2 and /3 2 . 
However, the FCS exposures could not be used to find r 2 
and /3 2 since no data was recorded in the first two seconds 
following the detector switch-on. This is also the reason 
why the hook response is not seen in Fig. 5. Furthermore, 
there is very probably a non-negligible timescale needed 
by the FCS to reach its final operating temperature. 
Therefore, the "fast" parameters (with timescales ^0.1 < 
r < ~ 5 s) were found using a self calibration technique 
(described in Sect. 4.2), operating on co-added repeated 
chopper sweeps crossing standard point source calibrators. 
During the self calibration the "slow" parameters were 
fixed to the values determined in the fits to the FCS ex- 
posures. Observations of standard point source calibrators 
(12 source- wavelength combinations for both the CI 00 
and C200 arrays) were analysed using the self-calibration 
technique. Each observation of a source provided several 
points for the determination of the "fast" detector con- 
stants, since individual subscans through the source could 
be processed independently. In all, self calibrations were 
performed on data spanning a range in detector illumi- 
nation of about 100. This allowed the dependence of t 2 
and /3 2 on illumination to be determined according to the 
parameterisation of Eqns. 8 & 9. Examples showing the il- 
lumination dependence of /3 2 and t 2 for the central pixel of 
the C100 array are given in Figs. 6c and 6d, respectively. 
The values of the 12 parameters determined for each pixel 
of the C100 and C200 detectors are given in Tables 1 and 
2, respectively. 
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Table 1. Model parameters for the CI 00 detector 
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4. The P32 algorithm 

In a standard reduction of ISOPHOT data using the in- 
teractive analysis procedures of PIA, the analysis is made 
in several irreversible steps, starting with input "edited 
raw data" at the full time resolution, and finishing with 
the final calibrated map. To correct for rcsponsivity drift 
effects, however, it is necessary to iterate between a sky 
map and the input data at full time resolution. This dif- 
ferent concept required a completely new data reduction 
package. In its development phase, and for the determina- 
tion of the detector model parameters, this data reduction 
package was run as a set of IDL scripts. Running in the 
PIA environment, P32Tools provides access to the main 
functionality of this data reduction package via a GUI in- 
terface (Lu et al. 2002). 

There are three basic elements to the concept of the 
P32 algorithm. The first step is signal conditioning, which 
aims to provide a stream of signal values, each with an 
attached time, from which artifacts such as glitches, dark 
current, and non-linearity effects have been removed, but 
which retains the full imprint of the transient response of 
the detector to the illumination history at the sharpest 
available time resolution. We will refer to such a stream 
of signals as the "signal timeline" . 

The second step is the transient correction itself, which 
is an iterative process to determine the most likely sky 
brightness distribution giving rise to the observed signals. 
This is a non-linear optimisation problem with the val- 
ues of sky brightness as variables. For large maps several 
hundred independent variables are involved. 

Lastly, there is the calculation of calibrated transient 
corrected maps. This entails the conversion from V/s to 
MJy/sr, flat fielding, and the coaddition of data from dif- 
ferent detectors pixels. Using P32Tools, this third stage 
requires the exporting of the transient corrected signal 
timeline into PIA at the SCP level, followed by the appli- 
cation of standard PIA processing. Here we just describe 
the new procedures for signal conditioning and transient 
correction. 



4.1. Signal Conditioning 

The processing steps for signal conditioning are as follows: 

1. Edited Raw Data is imported, after having been cor- 
rected in PIA for the integration ramp non-linearity. 

2. The specification of chopper phase in the data is 
checked and corrected where necessary. 

3. The pointing is established from the instantaneous 
satellite pointing information, unlike PIA which uses 
the commanded pointing. This also establishes the 
pointing directions on the slews between the spacecraft 
fine pointings. Internally defined "on target flags" are 
set for each spacecraft fine pointing direction, and the 
central direction and sampling intervals directions for 
the P32 "natural grid" (see Sect. 2.1) are established. 

4. The integration ramps are differentiated and the data 
is corrected for dark current and dependence of pho- 
tometry on reset interval, using standard calibration 
parameters taken from PIA. 

5. Systematic variations of the signal according to posi- 
tion in the integration ramps are removed from the 
data. This process is called "signal relativisation" . In 
contrast to the default pipeline or PIA analysis, this 
allows the use of all non-destructive reads, as well as 
providing a more complete signal timeline for use in the 
correction of the transient behaviour of the detector. 
The "signal relativisation" is applied such that signal 
variations due to the transient response behaviour are 
preserved. 

6. Random noise is determined by examining the statis- 
tics of the signals on each chopper plateau. This is nec- 
essary due to the high readout rates of the detectors in 
many P32 observations. Otherwise, the determination 
of noise from individual readout ramps, as done in a 
standard PIA analysis, would be insufficiently precise. 

7. Spikes are detected and removed from the data in 
a first-stage deglitching procedure. This operates at 
the full time resolution of the data given by the non- 
destructive read interval. 

8. A second stage deglitching procedure is performed, op- 
erating at the time resolution of the chopper plateau 
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Fig. 6. a) Values of the parameter f3\ for pixel 8 of the C100 array found from fits to the detector response to FCS 
exposures, plotted versus FCS illumination. The solid line is a fit of the function j3\ — (3io + flnS^ 2 with (3io — 0.96, 
/3n = —0.28, and f3\2 = 0.075. These are the parameters adopted for the model giving the default illumination 
dependence of /3i for this pixel, b) As a), but for the parameter t\ (in seconds). The solid line is a fit of the function 
T\ = Tio + TuS^J 12 with tio = 7.73s, r n = 11.60, and r 12 = —1.28. These are the parameters adopted for the model 
giving the default illumination dependence of t\ for this pixel, c) Values of the parameter /3 2 for the central pixel 
of the C100 array found from self calibration optimisation of the model on observations of standard point source 
calibrators, plotted versus illumination. The solid line is a fit of the function /3 2 — P20 + 02iS^ m with /3 2 o = 1.171, 
P21 = —0.870, and /3 22 = —0.0 145, where Sm um is the known illumination of the detector pixel due to the point 
source calibrator. These arc the parameters adopted for the model giving the default illumination dependence of 02 
for this pixel, d) As c), but for the parameter r 2 (in seconds). The solid line is a fit of the function t 2 = t 20 + ^l^mum 
with t 20 = 0.333 s, r 2 i = 0.381, and t 22 = 0.584. These are the parameters adopted for the model giving the default 
illumination dependence of t 2 for this pixel. 



(which comprise of at least 16 non-destructive read in- 
tervals) . This procedure (described by Pcschkc & Tuffs 
2002) removes any long lived glitch "tails" following 
the spikes detected in the first-stage deglitch. 

4.2. Transient Correction 

A data flow diagram summarising the iterative non-linear 
optimisation process to solve for the sky brightnesses sam- 
pled on the P32 natural grid is summarised in Fig. 7. The 
calculation to correct for the transient response behaviour 
is made separately for each detector pixel. The kernel of 
the procedure is to solve for the illumination correspond- 
ing to the signal on each chopper plateau. This is done in 



a constrained optimisation using a binary chop scheme for 
set values of the 12 detector parameters. The constraints 
set for the solution in illumination define the search range 
for illumination. They are set such that the solution for 
illumination should be positive, and less than 10 times the 
value of the highest uncorrected signal. In practice the il- 
luminations found on bright sources are never more than 
a factor of 5 higher than the uncorrected signal, so the 
determined photometry is not sensitive to the exact value 
of the upper limit in the search range. The values for S\ v 
and 52 P needed to solve for illumination on each chopper 
plateau are calculated from the response of the model to 
the preceding illumination history. After solving for the 
illumination on the current chopper plateau (in the repre- 
sentation in Fig. 7, to give the kth value Ik of the sequence 
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Fig. 7. Iterative scheme for transient correction (see Sect. 4.2) 



of illuminations in the "signal timeline" ) , the found value 
Ik is divided by the vignetting corresponding to the chop- 
per phase to obtain a new estimate of the sky brightess 
for the viewed direction on the P32 natural grid. This sky 
brightness estimate is then combined with any previous so- 
lutions for the same P32 natural grid pixel obtained with 
the same detector pixel from previous chopper plateaus, 
and the trial sky map (sampled on the P32 natural grid) 
is updated for this improved knowledge. The response of 
the detector pixel to the illumination history is then recal- 
culated and a solution is found for the illumination viewed 
on the next chopper plateau. This procedure is continued 
until all the signal timeline has been processed. Then a 
second pass through the data can be made (the "second 
major iteration" loop in Fig. 7), and so on, until the over- 
all solution for the sky brightness distribution no longer 
changes between successive major iterations. Any sky di- 



rections for which no solution can be found are flagged in 
an output mask, sampled on the P32 natural grid. 

The goodness of the solution is quantified through a 
value for % 2 , calculated from a comparison between the 
model fit and the signal timeline. The random uncertain- 
ties assigned to each signal sample in the timeline can 
optionally be recalculated from the self consistency of so- 
lutions for the illumination found in a given sky direction. 
Examples of model fits to the signal timeline are shown 
in Fig. 8 for the C100 detector and in Fig. 9 for the C200 
detector. 

The algorithm can also be used for datasets for which 
the readout rate was too rapid to allow all the data to be 
transmitted to the ground. In such cases the illumination 
history can still be determined at all times on fine point- 
ings, even where there are gaps in the signal timeline. On 
slews, the illumination history during the gaps are either 
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Fig. 8. Detail from a signal timeline from an observation of the standard calibrator Ceres, observed in the C105 filter. 
The green line shows the fit to the observed signal (shown in red) from pixel 1 of the C100 array, as found by the 
P32 algorithm for correction of the transient response behaviour of this pixel. 5 chopper sweeps, each comprising 13 
pointing directions ( "chopper plateaus" ) are shown, divided by vertical dotted lines separated in time by a duration 
of 6.1s for each sweep. The upper line in the figure shows the fine pointing flag; whereas the first chopper sweep shown 
was made while the spacecraft was slewing between fine pointings (flag value 0), the subsequent chopper sweeps were 
made on a single fine pointing (flag value 1). The remaining horizontal lines denote masks at the full time resolution 
of the data denoting glitches, readout status and whether or not a solution for the illumination could be found by the 
algorithm. The overall solution for the sky illumination for this detector pixel is given by the histogram in black. As 
described in the text, this overall solution is calculated from a combination of individual solutions for illumination on 
each successive chopper plateau, which are shown here as purple bars. For the very brightest sources, as shown in this 
example, the derived illuminations can be up to a factor of 6 brighter than the raw data for pixels in the C100 array. 
In this example the hook response is well seen in each chopper sweep at the 9th and 10th chopper plateaus. 



interpolated in time, or taken from the nearest previously 
solved direction on the P32 natural grid. 

For most observations using the C200 detector, and 
observations of moderate intensity sources made using 
the C100 detector, good fits to the entire unbroken sig- 
nal timeline can be obtained for each detector pixel. For 
observations of bright sources using the C100 detector, 
however, it is often necessary to break up the solutions 
into smaller segments, according to the distribution of 
the sources in the field mapped. This is best achieved by 
solving separately for the portions of the signal timeline 
passing through the bright sources. The most conserva- 
tive processing strategy, appropriate for observations us- 
ing the C100 detector of complex brightness distributions 
with high source to background contrast, is to break up the 
calculation into time intervals corresponding to individual 
spacecraft fine pointings. In this latter case, no processing 
of the slews is required. 



Several processing options are available for the inter- 
active processing, to optimise solutions for the particular 
characteristics of each dataset: 

Determination of detector starting state 

The default operation of the program is to assume 
that the detector is in equilibrium at the start of the 
observation, viewing the sky direction corresponding to 
the first detector plateau. This corresponds to the case 
Sip = Soi = Soo and Si p — 0. This will very rarely be 
the case, however, since an FCS measurement which may 
not be perfectly matched to the sky brightness will have 
been made just seconds prior to the start of the spacecraft 
raster. Therefore, the algorithm can optionally find a solu- 
tion for the detector starting state, as characterised by the 
values of Si p and S2 P immediately prior to the start of the 
first chopper plateau. This is done by solving for Si p and 
S2 P for a fixed S^ on the first plateau, where Soo is deter- 
mined from the last chopper plateau in the first spacecraft 
pointing to view the same sky direction. The timeline for 
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Fig. 9. Detail from a timeline of signal versus time from another observation of Ceres, this time observed in the C160 
filter using the C200 detector array. The key to the plot is as given in the caption to Fig. 8. This example shows the 
transition from a spacecraft fine pointing for which the chopper is sweeping through beam sidelobes to a fine pointing 
where the chopper samples the beam kernel. In this case the red line corresponding to the data is not seen during 
the intervening spacecraft slew, as it is exactly overplotted by the green line showing the fit to the data (see text). 
The chopper sweeps on the second fine pointing encompassing the beam kernel show a typical enhancement of the 
source-background contrast induced by the algorithm for correction of the transient response. 



the first spacecraft pointing is then repeatedly processed 
until a stable solution for Si p and S^p prior to the first 
plateau is found. This option was found to be particularly 
useful for data taken with the C200 detector. 

"Self calibration" of detector parameters 

In order to determine the primary detector parame- 
ters P2 and r 2 in the first instance, "a self calibration" 
processing technique was used by which any of the 12 de- 
tector parameters (Sect. 3.1) could be determined from 
observations of bright sources. The self calibration func- 
tionality has also been implemented in the publically avail- 
able P32Tools. This is because, even though the standard 
parameters of the detector model given in Tables 1 and 2 
generally give good results, some fine tuning of the param- 
eters may still yield improved fits to the signal timeline for 
some astronomical targets, particularly for targets with a 
high contrast to the background. 

Sources of any arbitrary brightness distribution and 
position can be used for the self calibration. However, in 
the initial determination of the detector parameters for the 
"fast" response of the detector, point source calibrators 
of known brightness were used, and the solutions for t 2 
and were additionally constrained by the requirement 



that the solutions for illumination were consistent with 
the known flux densities of the calibrators. 

The self calibration procedure works by simply repeat- 
ing the optimisation process for different subsets of the 12 
detector parameters, specified by the user. Each combina- 
tion of detector parameters produces an individual solu- 
tion for the sky brightness distribution, as well as a value 
for x 2 - -A- search is made in the parameter space of the 
detector parameters until a minimum is found in \ 2 ■ This 
is a lengthy process, however, making it advisable to per- 
form self calibrations on limited portions of the timeline. 
This is often chosen to correspond to a single spacecraft 
pointing where the chopper sweeps for the detector pixel 
being investigated pass through bright structure. The pro- 
cess can be further accelerated in cases that the repeated 
chopper sweeps through a source give the same repeated 
signal pattern. Then, the optimisation can be performed 
on an average of the coadded chopper sweeps performed 
on a single spacecraft pointing. This latter technique is 
referred to in P32Tools as a "composite self calibration" . 
Some examples of the useage of the self calibration proce- 
dure are given by Schulz et al. (2002a). 



Processing of slews 
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The calculation of the illumination history for data 
taken during the slews is less straightforward than for 
data taken during the spacecraft fine pointings. This is be- 
cause the instantaneous pointing directions during slews 
are effectively at random positions with respect to the 
P32 natural grid. Therefore, the use of solutions for the 
sky brightness distribution on the P32 natural grid found 
from analysis of previous fine pointings will in general lead 
to inaccuracies in the determination of the illumination 
history during the slews. There are two main approaches 
which can be adopted to counter this problem. 

The first is to solve for the illumination not over a 
chopper plateau, but for individual readouts in the signal 
timeline during the slew. This automatically provides a 
solution in which the model fit exactly matches the data, 
as there is only one fitted data point for each determined 
illumination. An example is given in the fit to the sig- 
nal timeline for the C200 detector shown in Fig 9. This 
approach is useful for observations of bright structured 
sources in which the variations in illumination with posi- 
tion on the slew on timescales of a chopper plateau exceed 
the signal to noise on individual data samples, ft is the de- 
fault processing option for the C200 detector. 

In cases of low signal to noise, and for all instances 
where not all data are transmitted to the ground, solving 
for the average illumination on each chopper plateau is 
a preferable technique. Otherwise, the values of <Si p and 
S2 P calculated just prior to the first chopper plateau on 
the next fine pointing can be wildly inaccurate, which can 
corrupt the solutions for the remainder of the signal time- 
line of the observation. Processing on a plateau by plateau 
basis is the default processing option to calculate the il- 
lumination history during slews for observations with the 
CI00 detector. Another reason for processing slews on a 
plateau-by-plateau basis for the C100 detector is that the 
procedure to solve for illumination for individual readouts 
is very susceptible to any residual glitches, and the C100 
detector is more prone to glitching than the C200 detector. 

5. Photometric performance 

In this section we give examples of the photometric per- 
formance achieved from maps made using the P32 algo- 
rithm. In all cases the standard calibration parameters 
used in the signal conditioning step prior to the correc- 
tion of the transient response behaviour and the subse- 
quent conversion from engineering units (V/s) to astro- 
nomical units (MJy/str and Jy) were taken from version 
8.1 of PIA. To illustrate the quantitative effect of the tran- 
sient correction we give results for datasets processed by 
the P32 algorithm both with and without the transient 
correction. In all cases, the corrections for the transient 
response behaviour were made using the standard detec- 
tor model parameters given in Tables 1 and 2. In Sect. 5.1 
we give examples of the effects of the transient correction 
on the derived spatially integrated flux densities. Where 
possible these integrated flux densities are also compared 
with flux densities from IRAS or from source models for 
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Fig. 10. The measured integrated and colour-corrected 
flux densities of the standard calibrators Ceres, Vesta and 
Neptune plotted against the nominal flux densities. The 
observations were done in various filters using the CI 00 
detector. The photometry derived with and without cor- 
rection for the transient response behaviour is shown with 
crosses and diamonds, respectively. The solid line shows 
the expected trend for equality between the measured and 
nominal flux densities. 

the standard calibrators. In Sect. 5.2 we give examples of 
the imaging performance. 

5.1. Integrated Flux Densities 

In general the corrections in integrated flux densities made 
by the algorithm depend on the source brightness, struc- 
ture, size, the source/background ratio, and the dwell time 
on each chopper plateau. The largest corrections are for 
bright point sources on faint backgrounds observed with 
the C100 detector. 

In Fig. 10 we plot the measured integrated and colour- 
corrected flux densities of standard ISOPHOT calibrators 
(the asteroids Ceres and Vesta, and the planet Neptune) 
versus the nominal flux densities of these sources (deter- 
mined as described by Schulz et al. 2002b). Since these 
same sources were used to derived the detector model, one 
expects a good statistical match between the nominal and 
the derived flux densities. Indeed, the integrated flux den- 
sities found after correction for the transient behaviour of 
the detector lie reasonably close to the solid line indicating 
an equality between the measured and nominal flux den- 
sities, though with some scatter (of around 15 percent). 
It can be seen that the flux densities derived from data 
uncorrected for the transient effects have been raised by 
factors ranging from ~ 2 for the fainter standards to ~ 5 
for the brighter ones. 

As an example of the photometric performance for 
fainter point sources we give results achieved for the faint 
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Fig. 11. The measured integrated flux densities of the 
faint standard star HR 1654 plotted against the predicted 
flux densities from a stellar model. The observations were 
done in various niters using the C100 detector. The pho- 
tometry derived with and without correction for the tran- 
sient response behaviour is shown with crosses and stars, 
respectively. The solid line shows the expected trend for 
equality between the measured flux densities and those 
predicted from the stellar model. 



standard star HR 1654. This source was not used in the de- 
termination of the detector model parameters, so it consti- 
tutes a test of the photometric performance of ISOPHOT 
in its P32 observing mode. The derived integrated flux 
densities, with and without correction for the transient re- 
sponse behaviour of the C100 detector, were compared in 
Fig. 11 with predicted flux densities from a stellar model. 
The corrected photometry is in reasonable agreement with 
the theoretical predictions. As expected, there is again a 
trend for observations with larger detector illuminations 
to have larger corrections in integrated photometry. 

A good linear correlation is also seen between in- 
tegrated flux densities of Virgo cluster galaxies (Tuffs 
et al. 2002a), derived from P32 ISOPHOT observations 
processed using the P32 algorithm, and flux densities 
from the IRAS survey (Fig. 12). The interpretation of 
these measurements constituted the first science appli- 
cation (Popescu et al. 2002) of the P32 algorithm. The 
ISOPHOT observations of Virgo cluster galaxies were fur- 
thermore used to derive the ratios of fluxes measured by 
ISO to those measured by IRAS. The ISO/IRAS ratios 
were found to be 0.95 and 0.82 at 60 and 100 /im, respec- 
tively, after scaling the ISOPHOT measurements onto the 
COBE-DIRBE flux scale (Tuffs et al. 2002a). 

The effect of the transient response behaviour on the 
integrated flux densities of very extended sources is illus- 
trated by the photometry of a few bright extended nearby 
galaxies. Thus, in Table 3, Col. 5, we give the integrated 
flux densities of NGC 6946, M 51 and M 101 taken from 
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Fig. 12. Integrated and colour-corrected flux densities of 
Virgo cluster galaxies measured in the ISO CI 00 filter 
versus the corresponding flux densities measured by IRAS 
in its 100 ^m band (taken from Fig. 7 of Tuffs et al. 2002a). 
The dotted line represents the relation ISO/IRAS=0.82. 
The star symbols are used for galaxies not fully covered 
by the spacecraft raster. 



Tuffs et al. (2003; in preparation). For NGC 6946 the in- 
tegrated flux densities are given here for the first time; 
those for M 51 and M 101 supercede the values given by 
Hippelein et al. (1996). The standard values for the de- 
tector parameters (as given in Tables 1 and 2) were used 
for the correction of the transient response. Table 3 also 
gives a comparison with the IRAS flux densities at 60 
and 100 /jm. In general, the ISOPHOT photometry yields 
somewhat lower integrated flux densities than the cor- 
responding IRAS photometry, similar to the ISO/IRAS 
ratios found for the Virgo cluster galaxies (Tuffs et al. 
2002a). As for the point sources, we also reduced the maps 
of these nearby galaxies without correcting for the tran- 
sient response of the detectors. We were thus able to de- 
rive the "correction factor" (Col. 6) indicating the ratio 
between the integrated flux densities with and without 
correcting for the transient. In comparison to the point 
sources of similar integrated flux densities, these correc- 
tions factors are moderate, varying between 1.07 and 1.45. 
This is because the timescale for the detector to traverse 
these extended sources becomes comparable to the time 
constants for the transient behaviour. The result therefore 
is a redistribution of flux density within the area of the 
galaxies along the scan direction, which goes some way to 
offset the loss of signal. We emphasise however that the 
remaining corrections are systematic, rather than random 
effects, and that the precise value of the correction factor 
for integrated flux density depends on the extent and mor- 
phology of the target. Furthermore, the correction factors 
for any point-like sources within extended sources, such as 
galactic nucleii or HII regions, will be much larger, com- 
parable to those found for the point source calibrators. 
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Table 3. Photometry of selected extended sources 



Source 


Size 


TDT 


Filter 


Integrated 


Correction 


Integrated 




(arcmin) 






Flux Density 
(ISOPHOT; Jy) 


Factor 


Flux Density 
(IRAS; Jy) a 


NGC6946 


9x7 


05302402 


C60 


111.5 ±6.5 


1.45 


146.5 


NGC6946 




05302401 


C200 


365.8 ±11.7 


1.11 




M51 6 


9x6 


17200618 


C60 


70.3 ±3.0 


1.30 


84.9 


M51 6 




17200618 


C100 


162.2 ± 12.3 


1.33 


192.6 


M51 6 




17200820 


C160 


180.9 ±23.9 


1.07 




M101 


21 x 10 


16500615 


C60 


61.1 ±3.2 


1.15 


92.9 


M101 




16500716 


C100 


204.5 ± 10.2 


1.12 


263.5 


M101 




16500716 


C160 


360.5 ± 49.0 


1.06 





a From Table 4 of Rice et al. (1988) 

b Together with companion galaxy NGC 5195 



5.2. Imaging Performance 

Fig. 13 shows an example of a brightness profile through 
HR 1654 at 100 /an, for data processed with and without 
the responsivity correction. Some 95% of the flux density 
has been recovered with the P32 algorithm. Without the 
correction, some 30% of the integrated emission is missing 
and the signal only reaches 50% of the peak illumination. 
The local minimum near 60 arcsec in the Y offset is a 
typical hook response artifact, where the algorithm has 
overshot the true solution. This happens for rapid chop- 
per sweeps passing through the beam kernel. This is a 
fundamental limitation of the detector model, which, as 
described in Sect. 3.1, does not correctly reproduce the 
hook response on timescales of up to a few seconds. This 
problem is particularly apparent for downwards illumina- 
tion steps. The only effective antidote is to mask the solu- 
tion immediately following a transition through a bright 
source peak. 
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Fig. 13. Brightness profiles (in MJy/str) along the Y 
spacecraft direction through the standard star HR 1654 
at 100 /im. The solid line represents the brightness profile 
obtained after processing with the P32 algorithm, while 
the dotted line shows the profile derived from identically 
processed data, but without correction for the transient 
response of the detector. 



Another effect of the inability to model the hook re- 
sponse is that the beam profile becomes somewhat dis- 
torted. This also has the consequence, that for observa- 
tions of bright sources, the measured FWHM can become 
narrower than for the true point spread function, as pre- 
dicted from the telescope optics and pixel footprint. We 
emphasise that, because of the beam distortion, the flux 
density of bright point sources cannot be found from the 
peak surface brightness, but instead should be determined 
by integrating the map. An example of an extremely bright 
point source reduced using the P32 algorithm is given in 
Fig. 14, depicting maps of Ceres in the C105 filter, made 
with and without the correction for the transient response 
of the CI 00 detector. 

Despite the limitations due to the lack of a proper 
modelling of the hook response, the algorithm can effec- 
tively correct for artifacts associated with the transient 
response on timescales from a few seconds to a few min- 
utes. This is illustrated in Fig. 15 by the maps of the 
interacting galaxy pair KPG 347 in the C200 filter (again, 
after and before correction for the transient response of 
the detector, respectively). The uncorrected map shows a 
spurious elongation in the direction of the spacecraft Y 
coordinate, which is almost completely absent in the cor- 
rected map. If uncorrected, such artifacts could lead to 
false conclusions about the brightness of FIR emission in 
the outer regions of resolved sources. Also visible are ar- 
tifacts at Y= ±310 arcsec. Similar artifacts are also seen 
at the same offset in other P32 observations made with 
the C200 detector, with position and amplitude relative 
to the central source independent of the filter used. These 
features are attributable to optical effects peculiar to the 
P32 mode, in which a small fraction of the light from a 
centrally positioned compact source is deflected away from 
the main beam into the vignetted part of the field of view 
for a certain combination of the offset in spacecraft fine 
pointing and chopper phase. Although these artifacts have 
nothing to do with the transient response behaviour of 
the detector, they are alleviated by the application of the 
P32 algorithm when solutions for detector illumination 
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Fig. 14. Top: Contour map of Ceres in the C105 filter 
after responsivity drift correction. The contour levels are 
5. 2" MJy/sr for all integers n with 1 < n < 10. The 
measured integrated flux density after background sub- 
traction is 98.4 Jy± 0.2 Jy (random) ±6% (systematic). 
The actual flux density of this standard calibrator (from 
a stellar model) is 109 Jy. As in Fig. 2, the the map pix- 
els of the P32 natural grid have been overlaid. Bottom: 
Contour map of Ceres in the C105 filter without correc- 
tion for the transient response behaviour of the detector. 
The contour levels are 5 . 2™ MJy/sr for all integers n with 
1 < n < 7. The The measured integrated flux density af- 
ter background subtraction is 28.0 Jy± 0.07 Jy (random) 
±4% (systematic) 

are converted into sky brightnesses, due to the vignetting 
correction made in this processing step. 

The image of M 101 in Fig. 16, made using the C100 
detector with the C100 filter, is given as a state of the art 
example of what can be achieved with a careful interactive 
processing of P32 data using the P32 algorithm. In addi- 
tion to the transient response corrections and a masking 
of residual hook response artifacts, a time dependent flat 
field has been applied. The spiral structure of the galaxy, 
with embedded HII region complexes and a component of 
diffuse interarm emission can clearly be seen. 
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Fig. 15. The interacting galaxy pair KPG 347 observed 
in the C160 filter. Top: after processing with the P32 al- 
gorithm. Bottom: with identical processing, except that 
the correction for the transient response behaviour of the 
detector has been omitted. The map pixels corresponding 
to the P32 natural grid (see Sect. 2.1) have been overlaid. 
In both maps, contours are logarithmic, at levels 1, 2, 4, 
8, 16, 32, 64 MJy/sr. The peak brightnesses are 126 and 
109 MJy/str for the maps with and without processing 
with the P32 algorithm, respectively. 
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Fig. 16. A 27 x 27 arcmin field containing the galaxy 
M101, as mapped in the C100 filter, after processing with 
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